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INTRODUCTION  AND  SUMMARY 


Under  the  previous  contract  AF  49 (638) -1507,  it  was  found  that  the 
heat  release  in  the  solid  phase  was  associated  with  the  driving  process 
of  unstable  combustion.  Since  most  of  the  propellants  of  interest  are 
based  on  ammonium  perchlorate  (which  has  been  shown  to  endow  propellants 
with  instability),  it  was  considered  necessary  to  identify  the  processes 
which  could  attenuate  acoustic  waves. 

To  date  our  study  of  loss  mechanism  has  shown  that  end  losses  are 
significant  and  enhance  both  axial  and  transverse  stability.  In  the  case 
of  axial  instability  the  enhanced  heat  transfer  to  the  end  closure  is 
being  examined  using  heat  transfer  gages. 

Our  theoretical  studies  are  continuing  in  an  attempt  to  integrate 
combustion  response  to  the  gas  dynamics  of  the  cavity. 
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EXPERIMENTAL  STUDIES 


Axial  and  Transverse  Mode  Dampening 

During  one  of  our  early  experiments  using  a  rectangular  slab  motor, 
a  core  became  malaligned,  and  as  a  result  the  center  rectangular  perfora¬ 
tion  was  significantly  out  of  alignment  with  the  axis  of  the  motor.  It 
was  decided  to  test  the  motor  to  ascertain  whether  this  was  a  significant 
variable  affecting  either  the  initiation  of  axial  mode  instability  or  its 
characteristic  response  in  the  motor.  On  firing  it  was  observed  that 
combustion  was  stable  even  though  pulsed  considerably  above  ihe  critical 
pressure  predicted  from  our  correlation  for  AP  propellants. 

It  was  postulated  that  stability  was  enhanced  because  the  traveling 
wave,  not  intersecting  normally  with  the  reflecting  surface  (head  end), 
lost  energy  after  reflection.  It  was  believed  possible  that  certain 
other  internal  design  factors  might  cause  similar  dampening.  One  pro¬ 
posed  design  was  the  introduction  of  a  head  end  which  was  deliberately 
out  of  alignment.  An  aluminum  wedge  5  inches  in  diameter  was  consequently 
attached  to  the  pulse  head  in  the  5-inch  x  50-ir.ch  motor  in  such  a  manner 
as  to  tilt  the  face  3CP  from  the  perpendicular  to  the  axis.  The  motor 
was  loaded  with  a  radial  burning  grain  of  PBAN  103  propellant,  which  has 
a  critical  threshold  pressure  of  600  psia.  During  the  firing  two  signifi¬ 
cant  observations  were  made.  First,  combustion  was  stable  to  980  psia 
before  the  third  pulse  drove  it  unstable;  during  instability  the  wave 
form  was  extremely  distorted  and  appeared  to  have  a  large  harmonic  con¬ 
tent.  Second,  in  all  of  our  past  testing  with  PBAN  103,  an  80/20  AP/ 
binder  propellant,  the  instability  invariably  changed  from  axial  to 
transverse  (tnis  normally  caused  higher  chamber  pressures  and  ruptured 
the  nozzle  safety  bolts).  In  contrast,  the  slant  nead  motor  remained 
in  the  axial  mode  throughout  burning. 

The  above  experiments  suggest  strongly  that  energy  losses  at  the 
ends,  in  particular  the  head  end,  contribute  significantly  to  stability 
in  both  the  axial  and  transverse  modes.  Heretofore  the  losses  at  the 
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ends  have  been  considered  insignificant  and  have  been  neglected  in  most 
of  the  analyses  of  over-all  stability  which  have  so  far  been  published 
(e.g.,  the  Hart-McClure  theory).  It  is  apparent  that  the  viscous  drag  on 
the  head  end  is  probably  a  significant  source  of  loss  for  transverse 
instability;  in  the  case  of  traveling  wave  instability,  our  data  (see 
next  section)  suggest  that,  in  relation  to  the  contact  time  of  the  shock 
front,  an  unexpectedly  high  rate  of  heat  transfer  occurs. 

The  influence  of  length  to  diameter  ratio  on  stability  is  obviously 
explained  by  the  fact  that  as  the  length  decreased  the  end  losses  became 
more  significant  relative  to  the  driving  from  the  burning  propellant. 

Heat  Transfer  During  Unstable  Operation  of  a  Motor 

Heat  transfer  rates  to  the  head  end  of  a  5-inch  x  40-inch  motor 
operating  in  both  the  stable  and  unstable  modes  have  been  measured. 

During  axial  mode  instability  the  heat  transfer  rate  was  approximately 
660  BTU  ft-2  sec-1 ;  during  stable  combustion  it  was  90  BTU  ft-2  sec-1 . 

This  high  rate  of  heat  transfer  at  the  head  end  suggests  that  the  inert 
head  end  serves  as  an  energy  sink  and  is  as  efficient  as  the  nozzle  in 
causing  appreciable  wave  pressure  losses.  At  this  time  it  is  postulated 
that  the  primary  mode  of  energy  loss  is  by  radiation.  The  thermal 
response  of  the  cool  head  end  is  low  compared  to  the  frequency  of  the 
instability,  and  therefore  the  end  is  seeing  an  on-coming  high  tempera¬ 
ture  shock  continuously.  This  same  loss  mechanism  must  operate  at  the 
nozzle  end;  however,  it  is  modified  by  several  factors. 

The  first  factor  to  be  considered  at  the  nozzle  end  is  the  influence 
of  the  nozzle  opening,  which  is  a  sink  for  mass  and  energy;  the  second  is 
the  two  different  surface  environments  interacting  with  the  wave.  In  our 
representative  test  mo+or  the  outer  annulus  behaves  somewhat  like  the 
head  end  except  for  the  shape  and  emissivity  (carbon  versus  steel  at  head 
end).  The  carbon  nozzle  would  nevertheless  present  a  rather  cool  reflec¬ 
tion  surface  and  absorption  surface  to  thermal  radiation  from  an  advancing 
shock.  The  nozzle  opening  area  is  '„ite  different  in  that  the  advancing 
pressure  wave  sees  a  high  temperature  shock  front  located  in  the  nozzle 
throat.  Radiative  heat  loss  in  this  case  would  be  expected  to  be 
minimized. 
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Wave  Growth  Studies  -  Influence  of  Aluminum  on  Gas  Dampening 


Several  experiments  were  conducted  in  5-inch  x  80-inch  slab  motors 
to  measure  the  wave  pressure  and  velocity  at  different  stations.  The 
propellant  variable  was  aluminum  content,  5%  versus  15%.  With  the  15% 

A1  loaded  propellant,  data  were  lost  due  to  Kistler  gage  instabilities 
induced  by  high  temperature;  neither  the  wave  pressure  nor  velocity  pro¬ 
file  could  be  deduced.  At  the  four  out  of  six  stations  which  functioned 
adequately  it  was  apparent  that  the  wave  pressure  was  significantly 
reduced;  this  test  has  to  be  repeated.  Unfortunately,  under  the  severe 
conditions  of  heat  transfer,  neither  motor  tube  nor  pressure  gages  last 
very  long. 

The  test  with  propellant  containing  5%  aluminum  was  more  successful 
in  that  the  wave  pressure  was  tracked  successfully.  However,  one  gage 
was  lost,  and  so  again  the  velocity  profile  could  not  be  deduced.  This 
propellant  generated  higher  wave  pressures  and  exhibited  higher  gain 
factors  than  did  the  80/20  AP/binder  propellant  reported  in  our  Annual 
Report,1  Contract  No.  AF  49 (638) -1507.  In  fact,  the  heating  produced  by 
this  instability  caused  a  case  burn-through  and  the  loss  of  two  motor 
cases.  We  plan  to  alleviate  this  difficulty  in  the  future  by  bonding 
two  strips  of  silica-loaded  buna-N  rubber  inside  the  casings  along  the 
slot , 

In  our  study  of  wave  pressure  and  velocity  in  the  5-inch-diameter 
motor,  it  was  reported  in  the  annual  report1  that  the  combustion  in  the 
5-inch  x  15-inch  motor  did  not  go  unstable,  and  we  were  unable  to  track 
the  pulse  beyond  one  cycle.  The  pulse  behaved  as  though  it  were  influ¬ 
enced  by  the  nozzle  attenuation  before  it  ever  reached  the  nozzle.  Thus 
it  was  postulated  that  there  might  be  a  difference  in  response  between  a 
hot  propellant  gas  and  a  cold  gas.  To  ascertain  whether  this  is  true, 
the  5-inch  x  15-inch  motor  was  cold-flowed  with  helium  and  pulsed. 


^■Capener,  E.L.  ,  R.  J.  Kier,  L.  A.  Dickinson,  and  G.  A.  Marxman,  Response 
of  a  Burning  Propellant  Surface  to  Erosive  Transients,  Final  Scientific 
Report,  A.F.  Office  of  Scientific  Research  Contract  No.  AF  49(638)-1507 
March  15,  1966 
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Helium  was  chosen  because  the  speed  of  sound  in  helium  and  propellant 
are  approximately  the  same,  3200  ft/sec.  The  helium  flow  rate  of 
0.23  lb/sec  was  sufficient  to  cause  choked  f^ow  at  the  nozzle  and  raise 
the  chamber  pressure  to  100  psig  using  the  same  nozzle  as  was  used  with 
propellant.  The  gas  was  introduced  at  the  head  end  of  the  motor  con¬ 
taining  an  inert  propellant  with  mechanical  properties  similar  to  those 
of  PBAH  103,  In  both  hot  propellant  gas  and  cold  helium  gas  the  pressure 
response  to  the  black-powder  pulse  was  nearly  identical.  In  neither  case 
could  the  wave  be  tracked  beyond  one  cycle. 

Spectral  Studies 

In  further  tests  with  our  2-inch  x  20-inch  slab  window  motor  we 
re-designed  the  pulse  tube  to  reduce  the  free  volume  at  the  head  end  and 
removed  the  entry  contour  from  the  nozzle  inlet  in  order  to  reduce  the 
attenuation.  (This  is  the  motor  originally  used  in  attempts  to  examine 
shock  wave  structure  photographically.)  It  is  planned  to  examine  spec¬ 
tral  transmission  and  emission  in  association  with  related  studies  in 
extinction. 
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THEORETICAL  STUDIES 


A  more  rigorous  derivation  of  the  stability  criteria  reported  under 
Contract  AF  49(638)-1507  has  been  performed.  The  conclusions  remain 
unchanged,  justifying  the  interpretation  that  heat  release  in  the  solid 
phase  is  not  conduci\e  to  stable  burning.  The  more  rigorous  derivation 
has  been  achieved  by  transforming  the  mathematical  description  of  our 
combustion  model,  which  includes  surface-coupled  reactions,  to  a  problem 
that  is  mathematically  identical  to  an  earlier  analysis  by  Denison  and 
Baum,2  which  excluded  surface-coupled  heat  release.  In  this  way  an  exact 
treatment  of  our  combustion  model,  based  on  the  Laplace  transform  tech¬ 
nique,  is  obtained  by  extending  an  available  treatment.  By  comparing  the 
results  with  those  obtained  previously  by  Denison  and  Baum  for  zero  sur¬ 
face  heat  release,  the  effect  of  surface-coupled  reactions  on  combustion 
instabij'ty  bounds  can  be  evaluated  directly.  Also,  by  comparing  our 
simple,  approximate  solution  with  the  exact  solution  (which  is  unwieldy 
and  can  be  presented  only  graphically  for  practical  purposes),  the  limi¬ 
tations  of  our  mathematical  simplifications  can  be  determined.  Ultimately 
this  may  lead  to  an  improved  approximate  treatment,  which  gives  results 
close  to  the  exact  solution  but  can  be  presented  in  terns  of  simple 
equations.  This  is  highly  desirable  in  formulating  a  practical  stability 
theory. 

By  comparing  our  model1  with  that  of  Denison  and  Baum,2  the  reader 
can  easily  verify  that  the  following  transformation  of  a  set  of  para¬ 
meters  defined  by  them  makes  the  two  models  mathematically  identical 
(the  subscript  DB  denotes  the  parameter  as  defined  by  Denison  and  Baum). 
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The  Self -Excited  Mode 

Equations  (1)  transform  Denison  and  Baum's  solution  into  a  Laplace 
transform  perturbation  analysis  of  our  combustion  model,1  in  which 
surface-coupled  reactions  are  permitted.  This  analysis  shows  that  the 
response  of  the  combustion  mechanism  to  pressure  perturbations  has  two 
separate  modes,  the  "self -excited"  mode,  and  the  steady  oscillatory  mode 
The  bounds  of  the  self -excited  mode  are  illustrated  in  Figs.  1-4  for 
selected  values  of  the  surface-coupled  exotherm  parameter,  0  ,  where 
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Note  that  Q  is  a  measure  of  the  fraction  of  the  total  heat  flux 
s 

into  the  grain  that  comes  from  surface-coupled  heterogenous  and  decom¬ 
position  reactions;  is  the  total  heat  flux  from  gas-phase  and  surface 

coupled  reactions;  Q  is  the  heat  flux  from  heterogeneous  (pressure- 

H 

sensitive)  reactions,  and  is  that  from  decomposition  (pressure-insen¬ 
sitive)  reactions.  When  9  =0,  all  energy  release  occurs  in  the  gas 

s 

phase,  and  the  corresponding  curves  in  Fig.  1  are  those  generated  by 
Denison  and  Baum.  For  typical  values  of  the  activation  energies  and 


surface  temperature,  E^/RT^  ~  20  and  E^/RT^ 


20.  Thus  if 
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Q  /Q™  ~  0.02  and  Q  /Q  ~  0.02,  so  that  surface-coupled  reactions  pro- 
HI  D  T 

vide  4  percent  of  the  total  heat  flux  to  the  grain,  the  values  of  9 

s 

is  about  0.4.  As  a  first  approximation,  one  can  conclude  that  the 
magnitude  of  9g  is  of  the  order  of  ten  times  the  fraction  of  total  heat 
flux  coming  from  surface  couple  reactions.  For  example,  Fig.  2  corre¬ 
sponds  to  a  propellant  in  which  about  5  percent  of  the  total  heat  release 
is  generated  in  the  surface  zone ,  while  close  to  1/3  of  the  total  is 
attributed  to  this  source  in  Fig.  4. 

For  each  value  of  9  the  self-excited  mode  is  separated  into  three 

s 

distinct  regions,  labeled  a,  b,  and  c  on  Figs.  1-4.  The  thermochemical 

constants  of  the  propellant,  which  appear  on  the  ordinate  and  abscissa 

of  Fig.  1,  determine  the  nature  of  the  response  of  that  propellant. 

In  region  a ,  the  pressure  perturbation  causes  the  self-excited  burning 

rate  mode  to  respond  as  an  exponentially  decaying  oscillation;  i.e., 

the  self-excited  mode  is  stable.  In  region  b,  the  response  is  an 

exponentially  increasing  oscillation,  and  in  region  c,  it  is  a  monotonic 

exponentially  increase  in  amplitude.  Thus,  for  a  given  Q  (or  percentage 

s 

of  surface-coupled  heat  release) ,  a  propellant  that  falls  in  either 

region  b  or  c  will  be  unstable.  From  a  different  viewpoint,  a  propellant 

with  given  thermochemical  constants  may  be  either  stable  or  unstable, 

depending  on  the  extent  of  the  surface-coupled  exotherm.  For  example, 

a  propellant  with  A  =  3,  1/ql  =  1.5,  would  be  stable  if  Q  =0.5,  but 

DB  DB  S 

unstable  if  Q  =  3 . 

s 

The  stability  limit  for  each  value  of  9  is  the  curve  that  separate 

s 

region  a  (the  stable  zone)  from  regions  b  and  c  in  each  of  the  figures. 

In  zones  b  and  c  a  pressure  disturbance  of  any  frequency  or  type  will 
induce  unstable  combustion  in  the  self-excited  mode,  with  an  oscillatory 
response  in  b,  and  a  monotonic  response  in  c.  In  region  a  the  amplitude 
of  the  self-excited  mode  of  response  tends  to  decay  with  time,  so  that 
this  mode  is  stable. 

In  Fig.  5,  the  stability  bounds  (limits  on  region  a)  for  9  =-0, 

s 

0.5,  1,  and  3  are  shown  together.  This  figure  illustrates  the  rather 
drastic  effect  on  stability  characteristics  which  arises  with  even  a 
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relatively  minor  appearance  of  exothermic  surface-coupled  reactions.  If 
as  little  as  5%  of  the  total  heat  release  occurs  in  such  reactions,  the 
extent  of  the  stable-operating  zone  is  greatly  diminished  relative  to 
that  with  gas-phase  reactions  alone. 

The  corresponding  stability  limit  derived  from  our  previous  approxi¬ 
mate  analysis1  is  also  shown  in  Fig.  5  for  each  of  the  indicated  values 
of  0g.  In  general,  the  approximate  treatment  agrees  with  the  exact  one 

for  small  values  of  A  ,  and  the  agreement  extends  to  larger  values  of 

DB 

A  aS  ®s  increases<  In  cases,  the  approximate  analysis  tends  to 
overestimate  the  zone  of  stable  operation. 

In  almost  all  practical  propellants  the  heat  release  in  surface- 

coupled  reactions  probably  represents  at  least  10  percent  of  the  total, 

so  that  0g  >  1.  Moreover,  for  most  propellants  the  decomposition 

activation  energy  is  such  that  A  <  5,  and  the  gas-phase  flame  tern- 

DB 

perature  and  activation  energy  are  such  that  l/<v  <  2.  Therefore, 

in  the  regime  of  greatest  practical  importance  the  approximate  analysis 
is  relatively  close  to  the  exact  one. 

The  Steady  Oscillatory  Mode 

In  addition  to  the  self -excited  mode,  which  either  grows  or  decays 
according  to  the  criteria  described  above,  there  is  a  steady  oscillatory 
mode  of  burning  rate  response  to  pressure  perturbations.  For  a  given 
value  of  0g,  this  mode  is  important  in  region  a  of  Fig.  1,  the  region 
in  which  the  self -excited  mode  decays  and  the  propellant  is  technically 
stable.  The  steady-oscillatory  mode  is  a  constant-amplitude  oscillation 
whose  amplitude  depends  on  the  perturbation  frequency  and  on  the  param- 
eters  ADfi  and  o^. 

This  behavior  is  illustrated  for  the  specific  case  0g  =  0.5  in 
Fig.  6.  Along  the  solid  curves  the  ratio  of  maximum  burning-rate  ampli¬ 
tude  to  perturbation  amplitude  is  constant  and  has  the  value  shown.  The 
perturbation  frequency  required  to  give  this  response  is  indicated  by 
the  dashed  curves;  any  other  perturbation  frequency  will  result  in  a 
lower-amplitude  response.  The  solid  curve  labeled  co  corresponds  to  the 
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limit  line  dividing  regions  a  and  b  in  Fig.  2.  The  zone  beyond  this 
curve  is  region  b,  in  which  a  pressure  disturbance  of  any  character¬ 
istic  frequency  will  induce  an  unstable  response  in  the  seif-excited 
mode  as  described  above. 

It  is  evident  from  Fig.  6  that  even  with  a  propellant  that  falls 
well  within  the  stable  zone  (region  a  of  Fig.  2  if  0^  =  0.5),  a  per¬ 
turbation  of  a  certain  frequency,  or  a  pressure  decay  whose  character¬ 
istic  time  coiresponds  to  that  frequency,  may  cause  an  undesirable  os¬ 
cillatory  response  of  greatly  magnified  amplitude  unless  the  motor  has 

a  high  damping  factor.  For  example,  a  propellant  with  0  =0.5,  A  =  3, 

S  DB 

and  l/«  =  1.5  is  technically  stable  according  to  Fig.  1.  However,  if  a 

DB 

pressure  disturbance  of  frequency  u)  =  2r2/K  is  imposed,  the  amplitude  of 

the  burning  rate  response  is  more  tnan  four  times  that  of  the  pressure 

disturbance,  assuming  the  order  of  the  gas-phase  reaction  n  2.  Thus, 

owing  to  this  characteristic  of  the  steady-oscillatory  mode  to  be  truly 

stable,  a  propellant  must  fall  well  inside  the  limit  line  between  regions 

a  and  b  corresponding  to  that  propellant’s  characteristic  value  of  0  . 

s 
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CONCLUSION 


Theoretical  and  experimental  studies  during  this  program  have  led 
to  the  inclusion  of  important,  previously  ignored  phenomena  in  a  com¬ 
bustion  model  that  is  reasonably  realistic  (consistent  with  available 
kinetics  information),  yet  mathematically  tractable.  The  early  exper¬ 
imental  studies  suggested  that  the  difference  in  solid-phase  decomposi¬ 
tion  reactions  for  different  composite  propellants  might  lead  to  greater 
fractional  heat  release  at  the  surface.  Through  an  approximate  analysis 
of  the  combustion  model,  it  has  been  shown  that  the  distribution  of  en¬ 
ergy  release  between  gas-phase  and  surface-coupled  reactions  may  play  a 
critical  role  in  determining  the  stability  bounds  for  a  solid  propellant.1 
The  stability  criterion  derived  from  the  approximate  analysis,  when  pre¬ 
sented  in  terms  of  pressure  and  burning  rate,  agrees  qualitatively  with 
experimentally  observed  stability  limits.1  A  direct  quantitative  com¬ 
parison  cannot  be  attempted  without  precise  kinetics  data,  which  are  not 
available,  of  course,  for  the  complex  propellant  combustion  mechanism. 

To  separate  implications  of  the  combustion  model  from  possible 
limitations  of  the  approximate  mathematical  analysis  of  that  model,1  the 
foregoing  comparison  with  a  more  exact  analysis  was  essential.  In  gen¬ 
eral,  this  comparison  has  confirmed  our  previous  conclusions1  qualita¬ 
tively,  and,  to  a  rather  surprising  degree,  quantitatively  as  well. 

The  analysis  summarized  above  indicates  that  when  even  a  very  small 
fraction  of  the  total  heat  release  occurs  in  surface-coupled  reactions, 
the  likelihood  of  unstable  combustion  is  greatly  enhanced.  The  limita¬ 
tions  indicated  by  the  present  analysis  m  this  respect  are  somewhat 
more  severe  than  those  shown  by  the  approximate  treatment.  The  more 
exact  analysis  also  shows  the  possibility  ox  a  high-amplitude  burning 
rate  response  with  certain  pressure-disturbance  frequencies,  even  when 
the  propellant  is  technically  stable  (in  the  sense  that  there  is  no 
"runaway"  burning  rate  response) .  This  phenomenon  has  the  effect  of 
broadening  the  limit  line  into  a  narrow  zone  that  separates  clearly 
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stable  propellants  from  those  that  always  tend  to  be  unstable.  The  re¬ 
sponse  of  the  propellants  that  fall  within  this  dividing  zone  depends 
critically  on  the  frequency  of  the  pressure  disturbance. 
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FUTURE  WORK 


It  should  be  noted  that  the  experimental  studies  are  determining 
the  magnitude  of  the  losses  which  should  be  incorporated  into  a  unified 
treatment  of  the  propellant  response  and  the  gas  dynamics  of  the  cavity. 
Additional  work  is  proceeding  using  adiabatic  self -heating  and  DTA  exper¬ 
iments  to  elucidate  the  different  characteristics  of  solid-phase  reaction 
for  these  propellants  which  exhibit  contrasting  stability  behavior. 

During  our  future  theoretical  studies  we  shall  endeavor  to  modify 
the  approximate  solution  previously  developed1  so  that  it  is  quantita¬ 
tively  closer  to  the  more  exact  treatment  outlined  above.  The  latter  is 
difficult  to  express  in  other  than  graphical  form,  and  is  therefore  in¬ 
convenient  for  the  purpose  of  obtaining  an  over-all  description  of  in¬ 
stability  bounds  in  rocket  motors.  On  the  other  hand,  it  should  be 
possible  to  combine  equations  for  the  transient  gasdynamicT  of  the  com¬ 
bustion  chamber  with  a  suitably  modified  approximate  analysis  of  t.^e 
combustion  model,  expressed  in  equation  form  and  encompas  mg  all  the 
important  features  mentioned  above,  to  accomplish  this  objective. 
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FIG.  3  STABILITY  BOUNDS  OF  SELF“E)'CITED 
MODE  FOR  0S  =  1 
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FIG.  4  STABILITY  BOUNDS  OF  SELF-EXCITED 
MODE  FOR  0S  =  3 
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